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We present a unified picture to explain the spectral states of BH binaries based on the two-flow 
model developed in our team since several years and already applied with success to AGNs. In 
our view, the central regions have a multi-flow configuration consisting in (1) an outer standard 
accretion disc (SAD) down to a transition radius rtr, (2) an inner magnetized jet emitting disc 
(JED) below rtr driving (3) a self-collimated non relativistic MHD jet surrounding, when condi- 
tions for pair creation are met, (4) a ultra relativistic pair beam. Large values of r,r correspond 
to Hard states while small values correspond to Soft states. In between these extremes, in the 
high intermediate state, r„- can reach values that switch on a pair cascade process giving birth to 
ultra-relativistic pair blobs that explain the superluminal events. In this model the accretion rate 
but also the disc magnetization play important roles, the latter being necessarily in the range 
0.1-1 for the MHD jet to exist. Then, with simple assumptions on jU, we propose an explanation 
for the hysteresis behavior observed in microquasars during their outburst . We also discuss the 
nature, in our framework, of the X-ray corona. While it is commonly believed to be at the base 
of the jet it cannot be in the jet itself, its optical depth being to small to be consistent with the 
X-ray spectra observed in the Hard state. We show that it is expected to form naturally in the 
JED itself. Finally we investigate the nature of the X-ray variability in the hard state by assuming 
an oscillating hot plasma whose temperature and density vary locally through the propagation 
of a magneto-sonic wave. The variable comptonized spectra are computed through Monte-Carlo 
simulation. Preliminary results are discussed. 
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1. The two-flow model 

The two-flow model consists in combining a relativistic flow along the axis of a non-relativistic 
or mildly relativistic MHD flow coming from the accretion disk (cf. Fig. [l]). When this model has 
been proposed in the eighties (mostly with a description of the interaction of the cold relativistic 
beam with an ambient medium, e.g. pO|]), the arguments were rather poor: the phenomenology of 
hot spots do not require a relativistic flow, whereas relativistic sporadic ejections are observed only 
at short distance from the source; the large scale jet does not seem to lose significant energy during 
its travel, and non-relativistic shock acceleration satisfactorily explains the hot spot synchrotron 
spectra. So the large scale non relativistic -jets and the short scale relativistic-jets could be generated 



differently. First as trophy sical applications of the model have been published in 1989 [144|, |38, 



But more arguments have been proposed later [g2[ g3|] in relation with X-ray and gamma ray 



emissions. The interest of the two-flow model relies on the following main arguments: 

1. The relativistic jet alone would radiatively cool rapidly and be relaxed by Compton drag 



[]42[]. The heating of the relativistic jet particles through MHD wave interaction with the 
non-relativistic flow suppresses or at least strongly weakens this effect. 
2. The relativistic jet does not self-collimate [^,37]. The non-relativistic MHD flow, which self- 



collimates, can confine it and manages an ecological niche for it, because of the centrifugal 
barrier inherently associated with the MHD acceleration of the non relativistic -jet. 
The extraction of rotation energy from a spinning Black Hole does not generate a power 
strong enough to account for the large power of FR2 jets. Jets powered by accretion are 



more powerful, but are mildly relativistic [37]. 



4. The relativistic-flow is hkely in the form of flaring clouds, rich in electron-positron pairs, 
channelled by the non relativistic-jet, and involving a smaller power. It probably does not 
travel beyond kpc distance from its source. 

A new version of the two-flow model has been developed since 1991, where the relativistic-jet 



is a hot pair plasma, energized by turbulence from the non relativistic-jet ||2j, ^ and experience a 



Compton Rocket effect. Satisfactory fits of the spectra of 3C273 and 3C279 has been obtained with 



a natural break at a few MeV [|29[]; the stratification, which appears hardly avoidable (e.g. []24|]) 
plays a major role in the analysis, which deeply differs from homogeneous or slab models . Later, 
a study of the micro-instabilities of such a relativistic-jet have been investigated [0O|]; these insta- 
bilities differ from the one investigated in the previous version, by the fact that they are triggered 
by the streaming of the cold ambient medium in the hot relativistic-jet, instead of considering the 
relativistic-jet as a cold beam exciting waves by pervading the ambient medium. This is an interest- 
ing circumstance of plasma turbulence excitation, because the streaming instability is driven by the 
radiation force and quasi-linear theory is not sufficient to describe the stationary state; non-linear 
mode couplings are necessarily taken into account, and, in particular, allow to estimate the heating 
of the pair plasma. 

2. Application of the two-flow model to microquasars 

The application of the two-flow model to microquasars is relatively natural. First, there is a 
large consensus on the fact that the physical processes occurring in AGNs and microquasars are 
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similar, the major difference being the global size of the system which depends linearly on the 
central black hole mass (from tens of solar masses in microquasars to millions or even billion 
of solar masses in AGNs). Then, it appears reasonable to apply a model that has been validated 
in the AGNs framework (cf. previous section) to the microquasars one. On the other hand, the 
importance of ejection processes in microquasars as well as the crucial links between accretion and 
ejections during their spectral state evolutions (deduced from accurate multi-wavelength studies, 



e.g. [Q,[13|]) imply that accretion and ejection processes have to be addressed all together in a global 
and consistent way. This is an important characteristic of our model where accretion and ejection 
are indeed treated in a self-consistent way (cf. below). The first self-similar solution of the complete 



set of equations of an accretion-ejection structure can be found in [ 17, 14]. Analytical computations 
and heavy numerical simulations have been developed since that time in different astrophysical 
contexts [jsj, ^, 15, 0, |l6|] and their validity have been confirmed by different authors (e.g. 



[26, 45|]). Finally, there are growing evidences of the presence of different jet velocities component 
in XrB, one moving at non or mildly relativistic speed while the other propagates in the ultra- 



relativistic regime |19, 2C, 12, 32]. While these results do not request univocally the presence of 



two different jet components, they are of course naturally explained by our model. 
2.1 General picture 

We assume that the central regions of BH XrB are composed of four distinct flows: two discs, 
one outer "standard" accretion disc (hereafter SAD) and one inner jet emitting disc (hereafter JED), 
and two jets, a non-relativistic, self-confined electron-proton MHD jet and, when adequate condi- 
tions for pair creation are met, a ultra-relativistic electron-positron beam. We also assume the 
presence of a hot thermal corona at the base of the jet that is more likely more likely part of the 
JED itself (cf. Sect. 4). A sketch of our model is shown in Fig. |] and the four dynamical compo- 
nents are discussed separately in more details in [|l8|]. This model provides a promising framework 
to explain the canonical spectral states of BH XrBs mainly by varying the transition radius rtr be- 
tween the SAD and the JED. This statement is not new and has already been proposed in the past 
by different authors (e.g. [11, ^, 27, 25]) but our model distinguishes itself from the others by 



the consistency of its disc-jet structure and by the introduction of a new physical component, the 
ultra-relativistic electron-positron beam, that appears during strong outbursts. 

The electron-proton MHD jet is collimated by large scale vertical field anchored somewhere 
in the accretion disc (the JED) and we assume that this large scale magnetic field has the same 
polarity everywhere. The presence of a large scale vertical field threading the disc is however 
not sufficient to drive super-Alfvenic jets. This field must be close to equipartition as shown by 



[]17p and [|14|]. The reason is twofold. On one hand, the magnetic field is vertically pinching the 
accretion disc so that a (quasi) vertical equilibrium is obtained only thanks to the gas and radiation 
pressure support. As a consequence, the field cannot be too strong. But on the other hand, the 
field must be strong enough to accelerate efficiently the plasma right at the disc surface (so that the 
slow-magnetosonic point is crossed smoothly). These two constraints can only be met with fields 
close to equipartition. 

An important local parameter is therefore the disc magnetization /i = / (HoPtot) where Ptot 
includes the plasma and radiation pressures. In our picture, a SAD is established down to a radius 
rtr where /i becomes of order unity. Inside this radius, a JED with ~ 1 is settled. At any given 
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Figure 1: Application of the two-flow model to microquasars. This is a sketch of our model showing its 
different dynamical components. A Standard Accretion Disc (SAD) is established down to a radius which 
marks the transition towards a low radiative Jet Emitting Disc (JED), settled down to the last stable orbit. 
The JED is driving a mildly relativistic, self-coUimated electron-proton jet which, when suitable conditions 
are met, is confining an inner ultra-relativistic electron-positron beam. The MHD power Pmhd flowing from 
the JED acts as a reservoir for (1) heating the inner pair beam (Pg+e-) ^nd (2) driving the compact jet (Pjet)- 
We also assume the presence of a thermal hot corona at the base of the jet, and which is more likely part of 
the JED itself (cf. Sect. 4). Field lines are drawn in black solid lines and the number density is shown in 
grayscale (logjQn/m^^). 



time, the exact value of r,^ depends on highly non-linear processes such as the interplay between 
the amount of new large scale magnetic field carried in by the accreting plasma (e.g.. coming from 
the secondary) and turbulent magnetic diffusivity redistributing the magnetic flux already present. 
These processes are far from being understood. For sake of simplicity, we will treat in this section 
rtr as a free parameter that may vary with time independently of the accretion rate (but see Section 3 
for a more realistic picture). In that respect, our view is very different from that of [11, 28] who 
considered only the dependency of the accretion rate M to explain the different spectral states of 
BH XrBs. 



2.2 Time evolution of BH XrBs 

The evolution with time of a BH XrB has been reported in Fig. ^. This is a synthetic Hardness- 
Intensity diagram (hereafter HID) as it is generally observed in XrBs in outbursts (e.g. [[T[ |l^]) . 
During such outbursts, the objects follow the A-B-C-D sequence shown on the figure before turn- 
ing back to A at the end of the outburst. We have over-plotted on Fig. ^ the different sketches 
of our model (this figure is clearly inspired by Fig. 7 of JTsI]) and we have plotted on Fig. ^ the 
corresponding SEDs expected in each spectral state. We detailed below the interpretation of the 
HID in our framework. 



• Ascending tlie Right Brancii: Let us start at a Low/Hard State (LHS) located at the bottom 
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Figure 2: Schematic Hardness- 
Intensity diagram as it is generally 
observed in XrBs in outbursts (this 
figure is clearly inspired by Fig. 7 
of []l3|]). During such outbursts, 
the objects follow the A-B-C-D se- 
quence before turning back to A at 
the end of the outburst. The jet line 
has been defined by Jl3| ] and sep- 
arates jet-dominated states to disc 
dominated ones. See Sect. 2.2 for 
more details. 



Hardness 



of the HID right branch (point A in Fig. Such state would correspond to a JED extending 
up to typically rtr ~ 10^ (in units of gravitational radii r^). This considerably lowers the 
emission from the inner radii of the SAD producing a weak UV/soft X-ray excess. The hard 
(1-20 keV) power-law component of photon index F ~ 1.7 is attributed to the warm thermal 
plasma at the base of the jet. The non relativistic MHD jet then produces the persistent IR 
and optically thick radio synchrotron emission (cf. the corresponding SED in Fig. ||). 

The Top Horizontal Branch: 

- Before the jet line : Arriving in B we assume that Vtr starts decreasing rapidly. Then, 
the JED undergoes an outside-in transition to a SAD. The BH XrBs enter the high 
intermediate state. The flux of the outer standard disc then increases while the JED is 
decreasing. Under such circumstances, the MHD Poynting flux released by the JED 
is still important (through the large M that characterizes this part of the HID) but the 
accretion-ejection structure (i.e. JED and MHD jet) itself fills a smaller volume, a direct 
consequence being a weaker emission of the thermal "corona" and the non-thermal 
MHD jet emission with respect to what it is while in the LHS. 

- At the jet line: The jet line has been defined by [|l^] and separates jet-dominated states 
to disc dominated ones. This correspond to the critical phase reaches by the system 
during its evolution along the top horizontal branch where the conditions for a strong 



pair production are fulfilled (cf. [18] for more details). In this case, we expect an ex- 
plosive behavior of the pairs, with the sudden ejection of blobs. The emission of these 
blobs, first in X and 7 ray and then, after a rapid expansion, in IR and radio, will proba- 
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Figure 3: The canonical spectral 
states of X-ray binaries . (I) Qui- 
escent state obtained with a low 
M and a large r,^: the Jet Emit- 
ting Disc (JED) occupies a large 
zone in the accretion disc. (II) 
Hard state with much larger M and 
smaller r,^: the pair creation thresh- 
old is not reached yet. (Ill) Lu- 
minous Intermediate state between 
the Hard and the Soft states: the 
high disc luminosity combined with 
the presence of a MHD jet allows 
efficient pair creation and relativis- 
tic acceleration along the axis, giv- 
ing birth to flares and superluminal 
ejection events. (IV) Soft state when 
there is no zone anymore within 
the disc where an equipartition field 
is present: no JED, hence neither 
MHD jet nor pair beam. 



UV/SoftX Hai-dX 



bly dominate the broad band spectrum, producing the hard X-ray tail and the optically 
thin radio emission present in this state (cf. the corresponding SED in Fig. The 
production of a series of blobs can even result in an apparently continuous spectrum, 
from radio to X/yrays. We note that, at the jet line, the rapid increase of the pair beam 
pressure in the inner region of the MHD jet, during e.g. a strong outburst, may dramat- 
ically perturb the MHD jet production and even makes it disappear in agreement with 
observations [13]. It is also worth noting that we do not expect all BH XrBs to reach 
a "pair production" phase along the top branch. Indeed, triggering the pair cascade 
process requires a large enough gamma-gamma opacity and hence a high accretion rate 
that could not be reached by all objects. 

- After the jet line: We assume that rt,- is still decreasing. We therefore expect the total 
disappearance of the JED and its MHD jet when rtr approaches the disc inner radius 
r,„, thereby also causing the end of the pair beam (if present). The inner regions of the 
BHXrB are a SAD with probably a strong magnetic activity. Indeed, the situation might 
be slightly more complex than a mere SAD because of the presence of a concentrated 
magnetic flux. No steady MHD ejection can be produced from the SAD but unsteady 
events could be still triggered. This is maybe the reason why this region in the HID 
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seems to harbor complex variability phenomena |[T], |36| ]. 

Descending the Left Branch: When XrBs reach the left vertical branch (point C in Fig. 
rtr is smaller than the inner disc radius r;„ i.e. the JED and the MHD jet have completely 
disappeared. The whole disc adopts therefore a radial structure akin to the standard disc 
model and we enter into the so-called soft state (also called thermal dominant state by 
where the spectra are dominated by strong disc emission (cf. the corresponding SED in Fig. 
|3|). The descent from C to D correspond to a decrease in intensity i.e. by a decrease of the 
accretion rate. This is the beginning of the fading phase of the outburst, r,^ keeping smaller 
than We note also that we still expect the presence of magnetic fields that may be the 
cause of particle acceleration responsible for the weak hard-energy tail generally observed 



in this state ([31, 46] and references therein). 



The Low Horizontal Branch: In D we assume that r,^ begins to increase again. Thus, ac- 
cording to this conjecture, there is an inside-out build up of a JED. Self-collimated electron- 
proton jets could be produced right away. This means an increase of rtr, the reappearance 
of the non-thermal MHD jet and the thermal corona and a decrease of the SAD emission. 
But, contrary to the Top Horizontal Branch, the accretion rate is now too low to allow the 
production of a pair beam. Consequently we do not expect superluminal motions during 
this phase. 



3. What can drive rtr ? 

From the previous sections, it is clear that the spectral shape of a BH XrB critically depends on 
the size of the JED relative to the SAD, namely r^^. As stated before, Tf,- is the radius where the disc 
magnetization jj. = Bj / (HoPtot) becomes of order unity. Thus, r,r depends on two quantities Ptot{r,t) 
and B^{r,t). The total pressure is directly proportional to M since Pfot = 

p,Q.k,h and M being respectively the disc mass density, the keplerian angular velocity, the disc 
scale height and the black hole mass. As a consequence, any variation in the outer SAD accretion 
rate implies also a change in the amplitude of the total pressure. But we have to assume something 
about the time evolution of the large scale magnetic field threading the disc. The processes govern- 
ing the amplitude and time scales of these adjustments of rtr to a change in M are far too complex 
to be addressed here. They depend on the nature of the magnetic diffusivity within the disc but 
also on the radial distribution of the vertical magnetic field. As already noted, up to now we have 
simply assumed that rtr and M are two independent parameters. However we can try to relate these 
two parameters with each other using simple but reasonable assumptions. 

If the disc magnetization /i has to be close to unity for an accretion-ejection structure to exist 
there is obviously some freedom on the acceptable value of this parameter. More precisely, analyt- 
ical and numerical computations show that }JL has to be in between ~0. 1 and ~ 1 and that it keeps 
roughly constant in the JED. Now, how is pL expected to vary in a SAD? Assuming that, in steady 
state, the poloidal field is mostly vertical i.e. no significant bending within the SAD, the radial 
distribution B^{r) in the SAD is provided by the induction equation: 

dB. 

Vm^'^UrB^ (3.1) 

or 
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where is the turbulent magnetic diffusivity. This equation has the obvious exact solution 

B, oc r"^'" (3.2) 

where = —rur/v,„ is the (effective) magnetic Reynolds number. In a turbulent disc one usually 
assumes that all anomalous transport coefficients are of the same order so that Vm — V,,, v,, being 
the turbulent viscosity. Since the (effective) Reynolds number = —rUrjVv = 3/2 in a SAD, one 
gets that any vertical magnetic field is naturally increasing towards the center. Now in a SAD of 
vertical scale height oc r^, the total pressure Ptot = scales as 

67rVv 

where Msad is the (constant) SAD accretion rate. Using Eq. (B^ we get \l oc r^^ with ^ = 



— 5 — 3/2. In a SAD ~ 3/2 and 5 is always close to unity (apart from the unstable ra- 
diation pressure dominated zone where 5 = 0). Of course, the real value of ^ critically depends 
on the magnetic Prandtl number i^,„ = Vy/Vm) but this result suggests that, in a SAD, one may 
reasonably expect /i to increase towards the center. In conclusion, a reasonable radial dependence 
of jJL should look like the one plotted on top of Fig. ^i.e. constant (but in the range 0.1-1) between 
r,„ and rtr in the JED and decreasing above rtr in the SAD. 



The hmited range of values of /i in a JED bring us to the following assumptions. First, at a 
given radius in the disc, we assume that: 

• a SAD-to-JED transition occurs when /x=l 

• a JED-to-S AD transition occurs when jx <0. 1 

On the other hand, jj. <x meaning that /i is controlled by M but also by the magnetic field 

strength. This ratio is not expected to be constant with M and we make the zero order approxima- 
tion that it anti-correlates with the accretion rate i.e. we assume that the variation of the accretion 
rate is the dominant cause of the variation of ^u. As we will see below, these different assumptions 
provide a simple way to explain the hysteresis behavior observed in microquasars. 

We have reported in Fig. ^ the light curve and hardness-intensity diagram of GX 339-4 during 
its 2002/2003 outburst as well as the expected evolution of the radial distribution of the disc 
magnetization jx. Starting from the initial quiescent state in A, where jj. is constant equal to 1 
in the JED and decreases in the SAD, we reach the highest flux level hard state B after a large 
increase of the accretion rate. Following our assumptions, from A to B, /i(r) simply decreases 
but keeps its global shape constant. Point B corresponds to the situation where /x(r„.) =0.1. A 
small increase of the accretion rate then produces a JED-to-SAD transition starting from the outer 
part of the JED, i.e. the transition radius rtr decreases. The system then reaches point C. This 
may correspond to two different situations. Either r,r reaches a critical values that switch on a 
pair cascade process giving birth to ultra-relativistic pair blobs, or r„. becomes to close to r,„ for 
the MHD jet structure to keep stable. Both cases lead to the same effect i.e. the disappearance of 



8 



A Unified paradigm for BHXrB 



P.-O. Petmcci 




50 100 150 200 250 300 350 400 
Time (MJD-52359) 





0.1 



SAD 



(D) 




•5n 



Figure 4: The expected evolution of the radial distribution of the disc magnetization jj. during an outburst 
where an hysteresis in the Hardness-Intensity diagram is observed. The light curve and HID are those 
followed by GX 339-4 during its 2002-2003 outburst See Sect. 3 for details. 



the MHD jet. In the former case, the pair blobs produce strong radio emission and cause possible 
superluminal events while in the latter case the transition should be less "explosive". In any case, 
the situation in C is expected to be quite unstable, the MHD jet structure oscillating between life 



or death while some pair blobs are randomly produced. As already said in Sect. this could 
explain the complex spectral and temporal behavior that characterize the upper left comer of the 
HID diagram. Anyway, from C to D we enter in a disc dominated state. Whatever the accretion rate 
does, and following our assumptions, the JED cannot re-appears unless becomes equal to I. 

In the case of GX 339-4, this corresponds to point D where M becomes sufficiently small for this 
condition to be fulfilled. From point D to A, a SAD-to-JED transition occurs starting from r,„ on 
and in A we turn back to the initial conditions, at the beginning of the outburst. Thus the hysteresis 
behavior is a direct consequence of the variation of /i during the outburst, variation controlled by 
the accretion rate. 

It is worth noting here that, due to the unknown dependence of the magnetic field with the 
accretion rate, SAD-to-JED and JED-to-SAD transitions are not expected to occur at the same 
accretion rate for different objects. Even in a given object the transitions can be different for 
different outbursts (like in the case of GX 339-4, cf. [Q]) since it depends certainly on the past 
history of the accreting system. 



9 



A Unified paradigm for BHXrB 



P.-O. Petmcci 



4. What about the X-ray corona? 

The radio/X-ray correlations found in XrB during their hard states (e.g. 10, ^ imply a 
strong link between the hard X-ray corona and the radio jet. It is now commonly believed that the 
corona is indeed the base of the jet and we show in this section that this corona is more likely a part 
of the JED itself. Indeed, as already said, self-consistent accretion-ejection simulations show that 
the magnetic field has to be close to equipartition with the total pressure in the JED i.e. va ~ Cg 
where va and c., are the Alfven and sound velocity respectively. In the jet however va is close to the 
speed of light (but still midly relativistic). The JED and the jet basis being threaded by the same 

/ hjED \ ^ 

magnetic field, this implies that the jet to disc density — — < where h jed is the JED 

njED \ r J 

hjED 

height scale. This is a relatively important constraint since simulations give at most e = ^0.1 

r 

and is more hkely < 0.01 . Thus the jump in density between the JED and the jet basis can be very 
large. Then it can be shown that the vertical optical depth of the jet basis is at most equal to: 

where m = Mc^ /LEdd and Hhase the height scale of the jet basis (in physical units). However fits of 
the X-ray spectra in the hard state require generally a corona optical depth of ~ 1 which implies too 
large jet basis size, completely inconsistent with observations (for example concerning the X-ray 
variability). 

Consequently, the hard X-ray "corona" is more likely in the JED itself. This is not unreason- 
able since we expect optically thin regions in the JED even at high accretion rate. This is a direct 
consequence of the torque produces by the jet in the accretion disc. First most of the accretion 
power goes in the jet so the protons in the disc are "cold" and the JED keeps geometrically thin (or 
at most slim). Secondly, the accretion velocity in the JED is much larger than in a SAD. Simula- 
tions show that m/^^ ~ u^^^ ^ m^'*^ then implying the vertical optical depth Xjed ^ ^sad- In 

njED 

conclusion, the radial profile of the vertical optical depth in the disc should look like Fig. g. For 
low accretion rate, the entire JED can be optically thin. For higher M, situations can be reached 




Figure 5: Top: The expected radial distribution 
of the vertical optical depth in the disc for a rela- 
tively low accretion rate. It strongly decreases 
inward at the transition between the SAD and 
the JED due to the strong decrease of the surface 
density in the JED (see Sect. 4). The JED is op- 
tically thin. Bottom: the same for a larger accre- 
tion rate. The middle part of the JED is optically 
thin while the inner and outer parts are optically 
thick. The optically thin region in the JED is ex- 
pected to play the role of the X-ray "corona". For 
even higher accretion rate, the entire JED is ex- 
pected to be optically thick. 
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where the JED is optically thin in its middle part but optically thick in its inner (close to r,„) and 
outer (close to r,^) parts. The last situation is very interesting since it could explain the presence 
of broad iron lines observed some times in hard state (e.g. [33, |^]) and requiring the presence of 
optically thick material close to the black hole. Work is in progress to better quantify these results. 



5. A toy model for the X-ray corona in the Hard state 

We have also investigated the origin of the high energy variability in microquasars focusing 
first on the Hard State variability. In our case (see previous section) this corona should look like 
a cylinder in the JED and we simply assume in this section that r,„ and r,^ are its inner and outer 
radii respectively. The out-coming spectrum then depends on the optical depth and temperature of 
this corona. The observed spectra in the hard state generally required t ~ 1 and kTg ~ 100 keV. We 
propose that the observed variability in this state is produce by the oscillation of the corona through 
external excitation in rtr producing compressive waves in the corona. This idea was already pro- 
posed by [|^] and we study more generally this problem through Monte-Carlo simulations. Here 
we focus on stationary compressive waves, the perturbations being totally reflected in r;„ . Perturba- 
tions dn{r,t) and dTg{r,t) of the corona density and temperature then produce spectral variability 
in the X-ray spectra. For low frequency oscillations, the density and temperature perturbation have 
the same sign in the corona. Thus we expect large effects on the output spectrum i.e. large rms. On 
the contrary, for high frequency oscillations, dn{r,t) and dTe{r,t) change sign between r,>, and rtr- 
In this case, multiple Compton scatterings smear the spectral variability effects on the out-coming 
spectrum. The larger the oscillation frequency the weaker the variability and the smaller the rms. 
Now, depending on the boundary conditions, resonances where dn{r,t) and dTg{r,t) reach large 
amplitude may exist producing QPOs in the Power Density Spectrum (PDS). We show in Fig. ^ 
the simulated Monte Carlo SED and PDS expected from a cylindrical corona with mean optical 
depth and temperature T = 1.4 and kT^ = 100 keV when this corona is excited by a external white 
noise and when r,„ and r,r are nodes of the compressive waves. The SED looks very similar to 
the common hard state SED. The PDS show a flat top profile at low frequency, decreasing at high 
frequency (above ~ {kTe/mpf-^ /rtr), in agreement with the PDS observed in this state. However 




10 100 0-0°l O-Ol 1 1° 

Energy (keV) Frequency (Hz) 



Figure 6: Simulated Monte Carlo SED and PDS expected from a cylindrical corona with mean optical depth 
and temperature i —\A and kTg — 100 keV when this corona is excited by a external white noise. In the 
PDS, we have added 3 different lorentzians to fit the broad band noise and the two major QPOs. 
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different QPOs corresponding to the different resonances of the system are also visible. While low 
frequency QPOs are known to be present in the hard state/High intermediate state, close to the 
cut-off frequency, the presence of multiple low frequency QPOs like in our simulations are rather 
uncommon (but see for example [43]). Works are in progress to include effects that could attenuate 
these QPOs (like wave transmission in r;„, Cabanac et al. in preparation). 
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